动态负载平衡的二维非结构网格DSMC并行算法研究

基于PC-CLUSTER群机并行体系结构与消息传递库MPI并行环境,研究了二维非结构网格DSMC并行算法。提出一类非结构网格动态分区策略,保证各子区域的分子数量大致相等,实现计算进程间的动态负载平衡。利用MPI库函数构造了符合DSMC并行原理的单步通讯法。采用单控制多数据流(SPMD)以及Master/Slave并行模式,设计了二维非结构网格DSMC整体并行算法。在程序的编制过程中,充分展现了Fortran90高级语言的主要特性,引人动态数组、指针、链表及派生类型数据。最后对过渡流域高超声速绕流进行了并行计算,数值试验的结果表明本文设计的并行算法可以取得较为理想的加速比。     

扩散抛物化(DP)NS方程组的意义及其在计算流体力学中的应用

本文首先讨论扩散抛物化（DP）NS方程组的早期研究工作：它的提出、数学性质、意义和在CFD中的应用，然后讨论扩散抛物化理论的一些新发展。这些新发展是对NS方程组数值计算进行物理分析的基础上得到的，其中包括NS方程组差分计算时，粘性剪切流对网格间距和格式精度的要求；粘性项只保留剪切项的广义扩散抛物化（GDP）NS方程组，它的性质和应用。由于高Re数流动在NS方程组的差分计算中，网格Re数彼此相差悬殊的特点，产生了计算离散单元守恒方程组的新的算法思路，即离散流体力学（DFD）算法。在DFD算法中需要同时计算三种不同的守恒方程组（Euler,DPNS和NS方程组）。本文讨论了DFD算法中需要同时计算三种不同的守恒方程组（Euler,DPNS和NS方程组）。本文讨论了DFD格式的构造、它的优点和应用。并以超声速绕前后台阶流动为算例，来说明GDPNS方程组的用处和DFD算法的优点。DPNS方程组、GDPNS方程组、DFD算法是高智提出的，对这些问题，他的合作者从理论、算法、数值验证和某些应用又作了系统的研究。     

先进复合材料格栅结构(AGS)应用与研究进展

 具有空间点阵的先进复合材料格栅结构（AGS）是一种新型高效的结构形式，在航空、航天结构中具有良好的应用前景。介绍了AGS优良的各项性能以及在美、俄等国航空、航天工程中的应用现状，着重描述了混合工艺法、模具膨胀工艺和拉挤 互锁格栅结构几种现行工艺方法，并对当前国内外针对格栅结构的力学性能和结构健康监测的研究进行了概述。最后提出了AGS进一步研究的工作展望。     

基于预测的计算网格负载平衡研究

计算网格中负载平衡有着重要的作用和显著的效果.预测节点性能是实现动态负载平衡的关键,本文首先引入预测网格节点性能的指数平滑算法,讨论了平滑系数的选取方法.然后设计一个负载平衡实验,实验证明使用基于该预测算法的动态负载平衡与静态负载平衡和不使用负载平衡时任务执行效率有明显的提高.     

无限插值方法生成三维网格

 本文研究了生成三维网格的无限插值理论,应用它构筑了绕翼身组合体的三维C-O型网格,这是一种代数方法,适用于不同拓扑结构的网格生成,具有程序编写简单、计算时间短的优点。在这种方法生成的网格上进行了二维与三维的数值实验,结果表明,它是一种有效的网格生成方法。     

多部件模型在全尺寸小型双函道涡扇发动机气流数值模拟中的应用

利用高阶单调Ｇｏｄｕｎｏｖ显／隐格式求解全尺寸小涡扇发动机Ｓ２流面上带粘性力项的非定常Ｅｕｌｅｒ方程组,获得了令人满意的稳定收敛解。对多部件模型方法在双函道小涡扇发动机气流数值模拟应用中存在的问题提出了建议。     

二维流热耦合数值模拟程序及方法研究

本文发展了一套二维流热耦合数值模拟程序,流体域中采用结构化网格,固体域中采用非结构化网格,流体域和固体域耦合面上的网格为非匹配连接。为了使该程序能够同时模拟低速不可压流动,程序中还加入了预处理方法。此外,本文还发展了一种能较好保持通量守恒的界面插值技术,来处理流体域和固体域耦合面上的非匹配网格。数值计算结果与实验值或商业软件结果的对比表明:本程序中的预处理方法有效地改善了低速不可压流动下的计算收敛性和收敛精度,使本程序能够同时模拟高速和低速不可压流动;扩压器和叶栅的流热耦合计算表明本文程序在流热耦合计算方面具有较好的模拟能力,同时也表明本文发展的界面插值技术的可靠性。     

子母弹抛壳过程非定常流动的数值模拟

利用弹簧近似和网格重构相结合的非结构动网格技术,耦合求解Euler方程及弹道方程,模拟飞行在46千米高空、0°攻角、马赫数6.5、壳片分离初速分别为6m/s,12m/s、壳片分离角速度分别为0°/s,50°/s,150°/s的子母弹抛壳动态过程。研究表明,抛壳初期,壳片与母弹周围形成复杂流动,在壳片内表面存在激波的生成、发展、减弱等非定常过程,导致壳片受到气动力随着分离距离增加表现为很强的非线性特征;壳片分离初始角速度对壳片姿态的影响不是线性的,尽管分离初始角速度增加到150°/s,由于转动惯量小、俯仰力矩大,壳片实际飞行姿态达不到设计要求,壳片攻角较小,气动阻力小,导致壳片伴随母弹较长时间飞行;增大分离初速是实现快速分离的有效措施。     

基于Euler方程的三维自适应笛卡尔网格应用研究

本文采用非结构自适应笛卡尔网格和有限体积法求解了三维Euler方程。在网格生成过程中，采用了基于模型几何数据、模型表面曲率和流场特征为基础的网格自适应生成及加密方法，采用模型特征恢复技术保证了计算模型与原始外形的一致。在计算中，将以中心差分为基础的Jameson有限体积法在三维笛卡尔网格中进行了推广。本文中对导弹、歼击机等复杂问题进行了数值实验，计算结果与风洞实验结果符合良好，并在工程实际中成功应用。     

双边界代数法在喷管内流场网格生成中的应用

用一种双边界代数方法生成二维网格，通过在超限插值函数中引入边界的一阶导数实现网格与边界是的正交，采用几何一致性条件保证物理区域不变形，发展一系列拉伸方程控制网格在边界或给定位置的间距。应用此方法生成了满意的喷管内流场网格，网格生成效率明显优于微分方程网格生成方法。